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Abstract 

A systematic numerical study of the classical solutions to the combined system consisting 
of the Georgi-Glashow model and the SO (3) gauged Skyrme model is presented. The 
gauging of the Skyrme system permits a lower bound on the energy, so that the solutions 
of the composite system can be topologically stable. The solutions feature some very 
interesting bifurcation patterns, and it is found that some branches of these solutions are 
unstable. 
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1 Introduction 



The physical motivation of the present work is to set the framework for a semiclassical ap- 
proach to describing the mechanism of monopole catalysis of Baryon number decay, proposed 
by Rubakov |I| and Callan 0. Here we are led by the work of Callan and Witten [Q, where 
the Baryon is described as the soliton of the Skyrme 0] model, in the backgroud of the U{1) 
Maxwell field of a monopole in the Dirac gauge. In the present work, the Baryon is again de- 
scribed by a Skyrmion, which in this case interacts with the full 50(3) non-Abelian Higgs model 
(the Georgi-Glashow model) so that the Skyrmion we consider is gauged with the full 5*0(3) 
group and interacts with the 't Hooft-Polyakov monopole We will refer to this model as 
Monopole-Skyrmion model (MSM). 

The most important difference between the present work and that of Ref. is that here we 
have two distinct topological charges - the first being the Baryon number of the Skyrmion and 
the second, the monopole charge. The energy of our composite system therefore has a topological 
lower bound consisting of a combination of these two charges of rather different geometric natures, 
the Baryon charge being a degree which cannot be expressed as a total divergence while the 
monople charge is a flux by virtue of being descended from the second Chern-Pontryagin class. 

The most important feature of describing the interaction of the Skyrmion with the gauge 
field in our Monopole-Skyrmion model is the presrciption of gauging the Skyrme field with the 
diagonal SU{2). This prescription was introduced in Ref. for the gauged Skyrme model, 
where no Higgs field and Higgs potential are present, and the resulting solutions were studied in 
Ref. 0, If. Most importantly, this gauging permits a lower bound on the energy of the gauged 



Skyrmion unlike in the case when the usual gauging is from the left, e.g. in Ref. |]TT| where 
there is no lower bound. In this paper, we shall refer to the models arising from the gauging 
prescription used in as gauged Skyrme models (GSM). Since the topological lower bounds for 
the GSMs were presented in detail in Refs. f^, ^, and because all that we need to know here is 
that these exist, we do not discuss them further here. 

Now the presence of a topological lower bound is not a sufficient condition for the existence of 
a topologically stable soliton. To illustrate this we refer to the graph of the energy versus Skyrme 
coupling constant in Fig. 1 for the GSM studied in Ij^, 0. Without being mathematically rigorous 
one can suppose that the branches Ags and i?gs correspond to solutions, which form local minima 
of the energy functional, whereas the connecting branch Ags correspond to solutions, which form 
sattlepoints. 

The bifurcation structure in Fig. 1 is quite different from that appearing in models where the 
Skyrme field is gauged as in [|lT| without a topological lower bound. The corresponding graph to 
Fig. 1 in that case is given in the work of Ref. |]T^ featuring only two branches as opposed to the 
three in Fig. 1. Of these two branches [jl2| only one corresponds to stable solutions, as expected 



from the work of Refs. |I3|. Thus a butterfiy pattern of bifurcations with two stable branches 
seems to be typical of of GSMs with topological lower bounds. We will find in our study of the 
Monopole-Skyrme model, that the butterfiy of Fig. 1 persits for some range of the parameters 
in the model. 

Having already presented the lower bound on the energy of Monopole-Skyrme model in 
Ref. PI, we do not repeat it here and proceed straight away to the study the solutions numerically, 
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Fig. 1 
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Figure 1: The energy of the gauged Skyrmion as a function of k in the region of the bifurcation. 



with a view to exposing some of their quahtative features that may be of some physical interest. 
The model we employ here is a slightly modified version of the one in Ref. |^. It features a 
particular interaction term between the Higgs field and the chiral field in such a way that it also 
fixes the asymptotic values of the chiral field subject to the finite energy condition, and assuring 
integral Baryon number. In Section 2 we present the model, discuss its spectrum, and discuss 
some some scaling properties which will become pivotal in the subsequent numerical analysis. In 
Section 3, we subject the system to spherical symmetry and give the classical equations to be 
integrated, the detailed numerical results of which we give in Section 4. In Section 5 we analyse 
the normal modes of the radial fluctuations around the solutions to the MSM, characterised by 
those values of the parameters for which the solutions display butterfly bifurcations, and verify 
that indeed the branches corresponding to A^s are unstable. A summary and discussions of our 
results are given in Section 6. 
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2 The Model 



The Lagrangian of the Monopole-Skyrme model is given by 



-^Tr {d^UU^D^'UW} - '^Tr { \d^UU\ D^U^] Id^'UW, D'^UW] } 

4'^ JgLL JL JJ 

+V{^,U))£r 



with 



l({$,f/} + {<|.,f/t} 



The field strength tensor of the 5^(2) gauge potential = A'^y is defined as 



d^A, - d,A^ + i[A^, A 



fll '-V\ 1 



(1) 

(2) 
(3) 



and the covariant derivatives for the Higgs field $ 
are defined as 



and the chiral matrix U = exp{z7r"ra} 



d^U + t[A^,U] 



(4) 
(5) 



respectively, g denotes the gauge coupling parameter, A the strength of the Higgs potential, i] 
the norm of the vacuum expectation value of the Higgs field, /,r the pion decay constant and k 
the Skyrme coupling parameter. The parameter (yf^$ characterises the direct coupling between 
the Higgs boson and the chiral matrix. 

The first three terms in (|l]) are the familiar Georgi-Glashow model which is characterised by 
the scale rj. The next two terms are the Skyrme model with covariant derivatives 0^11 allowing 
for an interaction of the chiral matrix with the gauge potential. The first of these terms introduces 
another scale, f-j^. The last term V{^, U) describes a direct coupling of the chiral matrix with 
the Higgs field, on which we will comment later. 

The Lagrangian (p is invariant under local SU{2) gauge transformations g. 



A 



u 



id,, 



The vacuum of the theory is given by the following constant configuration, 

A^ = , $ = r/rs , U = l . 
In order to identify the particle content of the model we expand around the vacuum. 



A, 



6al^ , <l> = rjTs + 6rra, U = {I 



(6) 



(7) 
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and insert into the equation of motion obtained from the Lagrangian (|I|). Neglecting quadratic 



terms in 5a°, 50° and 571°" we find 



a^a'^^vr" = glrfSn'' , a = 1, 2, 3 , 



d^d^'Sa'"'^ = , 
d^d'^Sa"^'' = g^ri^Sa"'" , a = 1,2 



where we have defined (5a|^ = 6a'^ H — 9^50^ea36, and used the gauge fixing conditions d^5a^^'^ = 

and dfj^Sd'^'"' = for a = 1,2. Thus we find for the masses of the gauge field, the chiral field and 
the Higgs field. 



mA = gr] , rn^ = giT<i-ri > ruh = 2y Xr] , (8) 

respectively. (We take all parameters g,gT,<s,,ri, X to be positive). Clearly, the mass of the chiral 
field stems form the interaction term of the chiral matrix with the Higgs field, which breaks the 
chiral symmetry. However, the spontaneous symmetry breaking does not lead to a mass splitting 
for the chiral fields. 

In the following we will motivate the potential term (j^) for the Higgs field and the chiral 
matrix. First consider the model without the potential term. /^From the finite energy conditions 
we find for the chiral field and the Higgs field in the limit r oo 

^ , D^U ^ . (9) 

Let us assume, that this condition is fulfilled for the radial component of the covariant derivatives 
and concentrate on the angular components. We decompose the Higgs field and the chiral matrix 
in the form 

$ = h(f) , = 0"r„ , with 0"0" = 1 , 

Uoo = cos/1 + isinfu , it = u"'Ta , with u"-u"- = 1 , (10) 

respectively, where h, 0" and /, are functions of the variables r, 6, (p. 

lYiom. the finite energy conditions follows that at infinity \h\ = rj, foo = const. However, 0" 
and if' may still be functions of the angular variables 9, y?. The conditions (P) now become 

£'a$oo = (5a0 + ^[^a,0]) =0, D^f/oo = slu/oo (9«U + m]) |^ = , (11) 



where a = 6,ip. The first condition can be fulfilled with A^\^ = (i(9„00 + ^^0) where 
A^J'^ is some function of 9, tp. The second condition can be fulfilled either with /oo = or with 



^aloo ~ (fc^aMM + A^^^iL^ , whcrc A^") is again some function of 6*, Lp. If the latter condition is 
fulfilled, /oo niay take arbitrary values. Hence, these configurations can be deformed continuously 
into configurations with trivial chiral matrix. 

In order to avoid this problem we have introduced the potential @) into the Lagrangian. 
Using the general decomposition ( PUf ) the potential can be written as V{^, U) = |^^^/^/i^ sin^/. 
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i. e. it couples the modulus of the Higgs field \h\ and the chiral function /, and does not depend 
on the "phases" 0", m". Consequently, the masses for the chiral fields, introduced by the potential 
at infinity, do not depend on the direction of the Higgs field in isospace. Furthermore, the finite 
energy condition now forces /oo = nvr, where n is an integer. 

The sum of the potentials in (|l]) has global minima at \h\ = rj, f = nn, and the matrix 
of the second variations has only non-negative eigenvalues at the global minima, provided the 
parameters A and are positive. 

In order to study the consequences of the two scales f.^ and rj we will take two points of view. 
First we will fix /^r and express all dimensionful quantities in units of /tt- Equivalently, we can 
fix rj and express all dimensionful quantities in units of rj. In each case, the ratio of the scales 
will enter the equations of motion parameter. 

We define the dimensionless quantities x = rfT^g , $ = $//^ , Vo = v/fn ■ Then the 
Hamiltonian becomes 



n[A,^,u] 



1 

--Tr {D^UU^D^UU^} Tr { \d^UU\ D^UU^] \d^UU\ D^UU^] ] 

4^ JgLL JL JJ 



^2 



Tr I ¥ 
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{$,[/} + {l>,f/t} 



(12) 



where we have defined X = \/ g"^ , k = kg , g.„q> = g-K<i>lg ■ In terms of masses Eq. (H) the 
parameters (^tt* and A can be expressed as (7^$ = m.„/mA and = 2mh/mA, respectively. The 
parameter 770 = f-K/v denotes the ratio of the scales. Apart form the last term the Hamiltonian 
( p!2D this is equivalent to the Hamiltonian studied before in 0|. The difference is, that we now 
consider r/o as a free parameter. In the limit 770 the minimum of the Higgs potential allows 
for a vanishing Higgs field. In this case we find the gauged Skyrme model considered before in 

Fixing the scale parameter rj we define x = rjgr, $ = rj^ and ^ = f^/rj. Then the Hamiltonian 
becomes 



n[A,^,U] = "1 J (^Tr{F^,F^''} + ^Tr{D^^D^'^ + ^Tr{{¥ 



C2 

-\Tr {d^UU^D^'UU^} - —Tr { \d^UU^, D^U^] Id^^UU^ D'^UU^] } 



4 

4 



Tr{^^ - 



-({|.,f/} + {<l,f/t} 



(Tx 



(13) 



with A = A/^f^ = A, K = = K and 5f,r<i. = S'tt^/s' = S'tt*- 

Comparing with the case of fixed scale /,r we find ^ = l/r^o, a; = r^ox, $ = $/?7o. Because the 
Hamiltonians (|12D and (0) are equivalent, we can obtain the properties of (p!2D from ([13|) and 

vice versa by using these relations. In particular, for the dimensionless energies E = and 



E 



9 
4nr] 



E 



Ti we have E = — . We will opt to work with (p!2[) in the following. 
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3 Static spherically symmetric equations 



The static spherically symmetric, purely magnetic Ansatz for the gauge field is P, 10 



Ao = 0,A, = c(x)|^, A, = (l-a(a:))^ + 6(x)|, = - sin ^ ((1 - a(a;)) | - , 

(14) 

where the su{2) matrices r^, a = r,6,ip are defined in terms of the Pauli matrices ti, r2, by 

Tr = sin6'(cosv9ri + siinpT2) + cosOr-^ , 
Tg = cos9{cosipTi + sin(pT2) — siu 6'r3 , 

Tip = — sin {pTi + cos ipT2 . 

The spherically symmetric Ansatz for the Higgs field is 

<! = h{x)Tr , (15) 

and for the chiral matrix 

U = cos/(x) +isin/(x)xr . (16) 
The Ansatz is form invariant under the U{1) gauge transformations 



g = exp^i-Trj , (17) 

where the gauge transformation function F is an arbitrary function of x. The gauge field functions 
transform as 

c{x) — > c{x) —xT'{x) , 

b{x) cos r(a;)6(x) — sin r(x)a(x) , 

a{x) cos r(2;)a(a;) + sin r(a;)6(a;) , (18) 

whereas the Higgs field function h{x) and the chiral function f{x) are invariant. To fix this gauge 
freedom we first impose the condition c{x) = 0, which still allows for global transformations with 
r = const. Further, we find that the functions a{x) and b{x) enter the Lagrangian only in the 
form a'^(x) + 6'^(x) and a^(x) + &^(x), which permits us to set b{x) = 0. 

With the Ansatz ([I^)-(p!6D restricted to the gauge fixing conditions and b{x) = the Hamil- 
tonian (|12|) becomes 

n[a,h,f] = ^ f (a'' + i^-f^ + ^ + a'h' + \{h'-v',rx' 



g J \ 2x' 

f I 2 ■ 2{- , A 2 2 ■ 2j!/ I-I2 , ^ f \ 

+ a sm / + 4k a sm /(/ H — j 



2 ' 2x2 



2 ^2 



^^l^h:'sin'f\dx . (19) 
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The differential equations for the functions a{x),h{x) and f{x) can now be obtained as the 
variational equations which extremize the Hamiltonian ([T9|), 



a" 

h" 
It 



ia' - I] 



x^ 



+ + sinV 



X^ 



f 



h' 

X 



-2- + /i 2- + 4A(/i2 - r]l) + gl^ sin^ 



8k a sin/ 



a cos/ 



0? sin^/ 



/'2)-2sin/a'/' 



+ 2 (sin/ cos/a^ — x/') 



+5f^$ sin/ cos//iV| 



x^ + 8/€^a^ sin^/ 



(20) 



These equations have to be solved due to boundary conditions which ensure regularity of the 
solution at the origin and finite energy, i. e. 



x = : a = l, /i = 0, f = n , 
X oo : a— >0, h ^ rjo , /— i>0. 



(21) 



4 Numerical results 

We have constructed numerically solutions of the model for several values of the parameters rjo, 
A, K and g.„q,. In particular we investigated the dependence of the solutions on the parameters 
and k. 

In Figs. 2a and 2b we show the energy as a function of rjQ for several values of k 

4vr/^ 

for A = and (^tt* = 0. For large values of r]o the energy is a monotonically increasing function 
of rjo- In the limit rjo oo the energy increases linearly with t]q, such that E/r]Q 1, i. e. the 
energy becomes equal to the Monopole energy (in units of Attt]). Indeed, this limit corresponds 
to the limit f^^O where the chiral field becomes trivial, U = —1 everywhere except at infinity. 

4.1 r/ <C A 

For small values of rjo the solutions develop bifurcations, corresponding to the 'butterfly' struc- 
tures in Fig. 2a, where for a fixed value of rjQ three solutions coexist. We observe form Fig. 2a 
that the bifurcations occur only for a finite range of the parameter k, k^J^ < k < where 
K^J), depend on A and g^,i.. For X = 0, g^.^. = we find ^ 0.374 and k^^^ ^ 0.4495. We 
demonstrate the details of the bifurcations in Fig. 2b for k = 0.4 as an example. This figure 
suggests that the branches B and A correspond to local minima of the energy functional, whereas 
the branch A corresponds to saddlepoint solutions. 

The bifurcation pattern looks similar to the bifurcation pattern found recently in the gauged 
Skyrme model [§. Indeed, the bifurcations in the Monopole-Skyrme model and the gauged 
Skyrme model are closely related to each other. In the limit i]q ^ the Higgs potential allows 
for a vanishing Higgs field. In this case we obtain the gauged Skyrmion model studied in Refs. 
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Fig. 2a 



Fig. 2b 
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Figure 2: (a) The dimensionless energy of the Monopole-Skyrmion solutions is shown as a 
function of rjo for several values of n for fixed A = and gn<i> = 0. (b) The same as (a) for fixed 
K = 0.4, A = and g-^is, = 0. r]l^^ and 770^'* indicate the values of rjo where the branches B and A 
merge with the saddlepoint branch A, respectively. 



1^, ^. Consequently, the Monopole-Skyrmion solutions should approach the gauged Skyrmion 
solution in the limit of vanishing rjo. In Fig. 3 we show in a 3D graph the energy of the Monopole- 
Skyrmion solutions as a function of rjo and k together with the energy of the gauged Skyrmion 
solutions. We observe, that indeed the energy of the Monopole-Skyrmion coincides with the 
energy of the gauged Skyrmion in the limit rjo ^ 0. 

In order to understand the behaviour of the solutions for small 770, let us first discuss the 
solutions of the Higgs-less gauged-Skyrme model. In this model the absence of the Higgs field 
allows two possibilities for the value of the gauge field function a{x) at infinity. For large values 
of K the function a{x) vanishes at infinity (branch -Bgs) whereas for small values of k the function 
a{r) approaches the value one at infinity (branches Ags and ^gs). These two cases correspond 
to the dashed and solid lines, respectively, plotted in the E — k, plane in Fig. 3. For a certain 
range of values < k, < three branches of solutions are present, reminding one of two local 
minima (i?gs and Ags) and a saddlepoint (Ags) of the energy functional. 

Now consider the limit ?7o — >■ for the Monopole-Skyrmion solutions. If /t < the Monopole- 
Skyrmion solution will approach the unique gauged Skyrmion solution represented by the solid 
line in Fig. 3, below k^. However, if < k < there are three different gauged Skyrmion 
solutions available. In this case each Monopole-Skyrmion solutions of the branches A, A and 
B approaches the corresponding gauged Skyrmion solution on the branches Ags, Ags and -Bgs, 
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Figure 3: The dimensionless energy of the Monopole-Skyrmion solutions is shown as a function 
of 770 for several values of k for fixed A = and g^^^ — together with the dimensionless energy 
of the gauged Skyrmion solutions as a function of and k% indicate the values of k where 

the branches B^s and Ags merge with the branch Ags, respectively. 



respectively. For k > the Monopole-Skyrmion solutions on the branches A and A Monopole- 
Skyrmion solutions cease to exist and the solutions on the remaining branch B tend uniquely to 
the gauged-Skyrmion (-Bgs) solutions represented by the dashed line above Kjin Fig. 3. 

The limit 770 is non-uniform for k < for the Monopole-Skyrmion solutions on the 
branches A and A. This is expected because the asymptotic values of the gauge field function 
a{x) for the Monopole-Skyrmion solutions and the gauged-Skyrmion solutions (branches Agg and 
Ags) are different. For the latter the function a{x) approaches the value one at infinity, whereas 
for the former ones a{x) vanishes at infinity. To illustrate the limit ?7o ~^ for small values of k 
we exhibit in Figs. 3a-3c as an example a sequence of field configurations of Monopole-Skyrmion 
solutions along the "butterfly" for k, — 0.4 in Figs. 2a and 2b. We follow the branch B down to 
r)Q^ (see Fig. 2b), continue with increasing 770 on the saddlepoint branch A up to 770^^ and finally 
we approach 770 = on the branch A. 
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Fig. 4c 
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Figure 4: (a) The gauge field function a{x) is shown for several values of 770 for fixed k — 0.4, 
A = and g-j^^ — 0. (b) The same as Fig. 4a for the scaled Higgs field function h{x)/r]o. (c) 
The same as Fig. 4a for the chiral function f{x). The different lines correspond to r)o — 1{B), 
r]o = 0.03(5), 770 = 0.03(i) , rjo = 0.03(A) and r]o = 0.0001 (^) from top to bottom. 



In Fig. 4a the profile of the gauge field function a{x) is shown for r^o = 1.0 and rjQ = 0.03 on 
branch B, for rjo = 0.03 on branch A, and, for rjQ = 0.03 and rjQ = 0.0001 on branch A. While 
for r)o — 1.0, 770 = 0.03 {B), and tjq — 0.03 (A), a{x) is a monotonically decreasing function of 
X, it develops a local maximum at some stage on the branch A while 770 increases. Passing to 
the branch A, now with decreasing 770- this local maximum persists. Along this path the local 
maximum a(a;inax) and its location and reach a(xinax = 00) = 1 as rjo tends to zero. 
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while the asymptotic region, where a{x) decays to zero, is shifted to infinity. Thus, in the hmit 
?7o — the gauge field function of the Monopole-Skyrmion solution tends to the corresponding 
function of the gauged Skyrmion solution for all x except at infinity. Therefore the convergence 
of the Monopole-Skyrmion on A and A to the gauged- Skyrmion solutions on Ags and Ags is 
non-uniform. 

In Fig. 4b we exhibit the profiles of the scaled Higgs function h{x)/riQ for the same parameters 
?7o like in Fig. 4a. Wc observe that h{x)/riQ is a monotonically increasing function of x for all 
770- However along the path of values of rjo described above, the magnitudes of the functions 
h{x)/r]o become progressively smaller on an increasing interval, while the asymptotic region, 
where h{x) approaches its vacuum value rjo, is shifted to increasing values of x. In the limit 
?7o ^ the function h{x) vanishes everywhere, signaling the merging of the Monopole-Skyrmions 
to the (Higgs-free) gauged- Skyrmions in this limit. 

In Fig. 4c we show the profiles of the chiral function f{x) for the same values of parameters 
rjo as in Figs. 4a,b. For all values of rjo the function f{x) in a monotonically decreasing function 
of x. In contrast to the gauge field function a{x) and the Higgs field function h{x) the chiral 
function f{x) docs not change considerably with tjq. 

The case discussed above in Figs. 4a,b,c pertains to k < k'^. For < k, < k'^^ there are 
gauged-Skyrmion solutions of both types, namely those on branches Ags and as well as on 
branch Sgs, so the convergence of the Monopole-Skyrmion solutions to the gauged-Skyrmion 
solutions can be both uniform and non-uniform . In that case the Monopole-Skyrmion solutions 
on the branches A and A approach the gauged-Skyrmion solutions on the corresponding branches 
A^Q and A^a, in a non-uniform way, for the same reasons as described above. However, the gauge 
field function of the Monopole-Skyrmion solutions of the branch B obey the same asymptotic 
behaviour as the corresponding function of the gauged Skyrmion solutions of the branch Sgs. 
For these solutions the convergence is uniform. 

For K > kJ, there is only one type of gauged-Skyrmion solution, namely those on branch 
5gS for which the asymptotic value of the gauge field function a{x) equals zero like for the 
Monopole-Skyrmion solution, and hence the convergence of these solutions as 770 —> is uniform. 



4.2 r/ > 

Let us now consider the case where the scale rj is much larger than the scale In Fig. 5 we show 
the field configurations for rjo = 2700 (solid lines) and for r]Q = 1 (dashed fines) for comparison. 
We observe that the gauge field function a{x) for rjo — 2700 approaches its asymptotic value at a 
very small distance from the origin. The same applies to the scaled Higgs field function h{x)/r)o, 
except for the long ranged tail, which is due to the power law decay for vanishing Higgs mass, 
i. e. for A = 0. The chiral function f{x), however, extends to larger distances from the origin. 
This is in contrast to the configuration for rjo = 1, where the change in the profile of all functions 
is roughly on the same interval. 

Note, that r]Q = 2700 corresponds to the case where the parameters /,r and are of the 
magnitude of the pion decay constant in low energy QCD and the vacimm expectation value of 
the Higgs field in the Weinberg-Salam model, respectively. For the energy of this solution we 
find AttE — 2700.013 47r, i. e. roughly the energy of the BPS Monopole (ATrrj). An appealing 
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Figure 5: The profiles of tlie gauge field function a{x), the scaled Higgs field function h{x)/rjQ 
and tlie scaled cliiral function f{x)/7i are shown for rjo = 2700 (solid lines) and r/o = 1 (dashed 
lines) for fixed n = 0.4, A = and gfTr* = 0. 



physical interpretation of this solution seems to be that a Monopole resides at the center of a 
baryon and dominates its mass. 

This can also be seen in a different way using gauge invariant quantities like topological 
charges. We define the topological Monopole charge density as 

PMP = ^Tr{F,,Dk^}e'^'' , (22) 

and according to Ref. |^ the gauge invariant baryonic charge density as 

1 

where we have defined for U = exp(i7r"r' 



e = cos Ivrl , r = sin |7r|— , = d.C + e^^^^A^C , A^' = , (24) 

|7r| 

with a,b,c,d run form 1 to 4 and a, /3, 7 run from 1 to 3. For the spherically symmetric Ansatz 



14[jlq) and with the dimensionless coordinate x the scaled charge densities become 



[hjl-a')]' , . [/+(l-2a^)sin/cos/r 
P^'= 4vrx^,o = • ^''^ 
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For 770 = 0.0001 and rjQ = 2700, with fixed k = 0.4, A = 0, g-^i^, = 0, we show in Fig. 6 the 
functions pmp = 4:Ttx^pmp (sohd fines) and pb = 4:'kx'^pb (dasfied fines) normafized by their 
respective maxima. The values of the normafization constants are (47ra;^pMp)max = 2.113 and 




Figure 6: Tfie profiles of tfie normalized topological functions pMp{x)/meix{pMp) (solid lines) 
and pB{x)/m&x{pB) (dashed lines) are shown for 770 = 2700 and rjo — 0.0001 for fixed k, — 0.4, 
A = and Qt^^ — 0. 



(47ra;2ps)max = 11080 for rjo = 0.0001, and (47ra;^pMp)max = 0.31 and (47rx^pB)inax = 0.113 for 
770 = 2700, respectively. 

/^From Fig. 6 we observe that for 770 = 2700 the location of the maximum of the function 
Pmp resides at a considerably smaller distance from the origin than the the location x^^^ of the 
maximum of the function pb- This confirms the interpretation as a Monopole inside a baryon. 
Note, that the function pb possesses an additional local maximum at a larger distance from the 
origin than the global maximum and a minimum with vanishing magnitude between the both 
maxima. We found that half the baryonic charge stems from the area below the first peak and 
half from the area below the second peak. This can be understood as follows. For distances larger 
than or near the location of the minimum Xmin the gauge field function a{x) almost vanishes. 
Setting a{x) — in Pb we find that the minimum corresponds to /(xmin) = 7r/2. Splitting the 
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integral of the baryonic charge density into two parts, 

^ — psdx + / pbcIx (26) 

we find for both parts the value 1/2. 

Let us next discuss the charge densities for rjQ = 0.0001. In this case we observe from Fig. 6 
that the peaks of pmp and ps are roughly at the same location. However, at a large distance 
from these peaks the function pmp possesses a second local maximum. We found, that the 
magnetic charge stems mainly from this second local maximum, whereas the contribution from 
the global maximum is marginal. This behaviour becomes plausible if we consider the profile 
of the gauge field function a{x) in Fig. 4a. For small values of x the function a{x) is close to 
one. Consequently, the magnetic charge density is small. For larger values of x the function 
a{x) develops a local minimum, this leads to the first peak of Pmp- When a{x) approaches 
again the value one, the function pmp again becomes very small. At large distances the function 
a{x) decays to its asymptotic value. This leads to the second maximum of Pmp- However, the 
magnetic charge density also depends on the Higgs field function h{x), and one would expect 
that the charge density has to be small if h{x) is small, i. e. in the region where the gauge field 
function possesses its local minimum, see Figs. 4a and 4b. This is indeed the case. Note, that for 
r)Q = 0.0001 the normalization constants of the function ps is several orders of magnitude larger 
then the normalization constant of the function Pmp- Thus, the Monopole charge density at the 
first peak is indeed very small compared to the baryon charge density. 

In analogy to the interpretation of a Monopole inside a baryon for large values of tjq, one 
could interpret the solutions for small values of ?7o as a baryon inside a Monopole. However, the 
gauge field is nontrivial at the location of the baryon and the Higgs field is not in the vacuum 
in this region. Thus the baryon would reside in a nearly symmetric phase. One may speculate, 
that this scenario might be interesting in respect to the decay of the baryon. 

4.3 Qt^^ > 

In most calculations we fixed the parameter (7^$ = 0. In view of the discussion in section 2 
this needs some clarification. The interaction term of the chiral matrix with the Higgs field was 
introduced into the model basically for ideological reasons as it fixes the value of chiral function 
at infinity, foo — (say). Then we assumed that in the limit Qt^^ — > the asymptotic value of 
the chiral function is stil fixed, and that the field configurations behave smoothly. Indeed, we 
found from our numerical analysis that this is the case and that the assumption is justified. 

Let us now discuss the case where (^Tr* is finite. As long as (^Tr* is small, the dependence of the 
solutions on the parameters r/o and k, does not change considerably. In particular the bifurcation 
pattern for small values of rjo as shown in Figs. 3 persists for small g^Tr*- The reason is simply 
that (7^$ enters the differential equations only as a factor of the Higgs field function h{x). In the 
limit rjo the Higgs field function vanishes and consequently the equations do not depend 
on Qt^^ in this limit. To discuss the more general case of finite rjo let us assume that for some 
parameters k, A, rjo solutions on the branches A, B and A coexist and form a butterfly in the 
E — rjo diagram for = 0. Then the butterfly will persist for small values of g.,^^. As g^^^ 
increases, the butterfly shrinks in size and disappears at a critical value of e. g. g^^ ~ 0.4 
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for fixed k = 0.4 and A = 0. Thus, the bifurcation pattern is similar to Fig. 3, if we replace k by 
^0) ''70 by g-ir^ and interchange the role of A and B. 

We now consider the case where (7^$ becomes very large at fixed parameters k, X,riQ. In 
the limit (^vr^ 00 the potential (0) becomes a constraint, which for the spherically symmetric 
Ansatz (|15D , (|16D becomes 

/i^(a;) sin^/(x) — >• as (^tt* 00 . (27) 

Note, that this constraint can neither be solved by a vanishing Higgs field function, h{x) = 0, 
because this violates the boundary condition h{x —>■ 00) = tiq, nor by a constant chiral function, 
f(x) = or TT, because this violates boundary conditions at the origin or at infinity. However, 
there is a third possibility. If the Higgs field function vanishes on the interval [0, Xq] and the chiral 
function vanishes on the interval [xq, 00] then the function h'^{x) sin^/(x) vanishes everywhere. 




Figure 7: The profiles of the gauge field function a(x), the scaled Higgs field function h{x)/rjQ 
and the scaled chiral function f{x)/ii are shown for (7^$ = 10*^ (solid lines) and g.„i^ = (dashed 
lines) for fixed k = 0.4, A = and t]q = 0.03 an the branch A. 



In Fig. 7 we show the field configurations for (7^$ = 10'^ and for g^^^ = for comparison for 
fixed parameters k = 0.4, A = and rjo = 0.03 on the branch A. We observe that the Higgs 
field function h{x) is indeed almost zero on the interval [0,a;o], with xq ~ 23.65, while the chiral 
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function /(x) changes continuously from /(O) = vr to /(xq) ~ 0. On the interval [xq, oo] the 
chiral function is almost zero and the Higgs field function changes continuously form h{xo) ~ 
to its asymptotic value h{oo) = rjQ. The figure suggests that in the limit (7^$ — * 00 the derivatives 
of the functions h{x) and f{x) will be finite but non-continuous at xq, whereas the gauge field 
function remains twice differentiable at xq. 

Taking into account the behaviour of the functions for large QT^q,, we observe from the differ- 
ential equation for the chiral function f{x) (|20|) , that the 7r$ interaction term will be almost zero 
for X < Xq. Assuming that h{x) increases linearly at Xq, we find for large g^,^ that for x > Xq the 
chiral function decays exponentially with the exponent ~ — S'tt*. Hence, we find the following 
scenario. For x < xq the interaction of the Higgs field with the chiral field vanishes, whereas for 
X > Xq the chiral field becomes increasingly massive. Consequently, the baryon is trapped inside 
the Monopole. On the other hand, because the magnetic charge density is proportional to the 
Higgs field, there will be (almost) vanishing Monopole density for x < Xq for large g.^is, and the 
Monopole is expelled from the baryon. 



5 Normal modes 

To show the instability of a solution of Eqs. (|20|) we determine the eigenvalues of the fluctuation 
matrix around that solution. The existence of a negative eigenvalue of a normalizable fluctuation 
mode indicates that a deformation of the solutions in the direction of this mode lowers its energy. 
Hence this solution can not be stable. Therefore, to show the instability of a solution, it is 
sufficient to find a normalizable fiuctuation mode with negative eigenvalue. For the discussion 



of the normal modes we will adopt the methods discussed in Refs. |T5l . 

Here we consider only radial fiuctuations around the solutions. We use the dimensionless 
coordinate x from the beginning. We introduce into the Ansatz small space-time dependent 
fiuctuations \l/a(a;)e*'^*, a = a,b, c, h, /, 

a{x) a{x) + ^„(x)e*'^* , 

b{x) b{x) + ^b(x)e^'"* , 
c(x) ^ c(x) + V2^e(a;)e^'^* , 

h{x) h{x) + ^^e'-* , 

X 

fix) ^ /(x) + ^^e-S (28) 

X 

and expand the Hamiltonian and Lagrangian to second order in the fiuctuation function \1/q,. As 
we are interested in the fiuctuations around the static solutions, we set 6 = 0, c = and extremize 
the Lagrangian in the background of the functions a, h, /, i. e. whenever second derivatives of 
these functions appear, they are replaced by the right hand side of the differential equations (pOf ). 
The form invariance of the Ansatz (0)-(H), expressed by Eqs. i\l7t }, ([TsD , reflects itself in 



the existence of normal modes with vanishing energy eigenvalue. These gauge zero modes obey 
the conditions 



= , = ^ (—) , ^^, = , ^f = . (29) 



72 V 



a 



17 



Because we are interested in the non-zero modes, we want to exclude the gauge zero modes. This 
can be done by imposing the conditions that the normal modes have to be orthogonal to the 
gauge zero modes with respect to the metric 

< ^, ^ >bc= J (^fe^fe + ^c^c) + . . . dx . (30) 

(We will give the complete form of the metric later). This leads to the condition on the functions 

K[^f„^,] = {x^,y -V2a^i, = . (31) 

To exclude the gauge zero modes we add ^K"^ to the Lagrangian, where /i is a Lagrange multiplier. 

/^From the system of differential equations we find, that the functions and couple to 
each other, but not to the functions "^a, ^/i and "iff and vice versa. Thus, we have two decoupled 
systems of differential equations, which can be solved separately. 



5.1 The system {^5*c} 

Let us first address the system {\E'b\E'c}- The differential equations become 



1 + _ „ /:zxa' — a 



x^ 



x^ 



^6, 



and the corresponding static Hamiltonian is 



9 

+2^2 



x^ 



l+2a%2 T/2 o^c^c 

—^l + + 2^^ 



(32) 
(33) 



x" 



X 



X 



- vl/,vl/^-^(vl/,vl/J 



X 



dx 



with 



'3a2 



bb 



fi' + sin'/ + 4k" sin"/ I /'^ + 



x^ 



9 • 1 r'^ 

a sm J 



x^ 



(34) 



(35) 



where we have set the Lagrange multipher equal to |. The boundary conditions for the 
functions and are given by 



x = : ^, = 0,^6 = 0, 



(36) 



For solutions of (|32|) , ( p3D we can evaluate the energy integral (|3^) by integration by parts 



and using 



and (^). We find for the dimensionless energy 

Ebc = ^'^ , 



47r/, 



-n 



bci 



(37) 



if we assume that the functions are normalised with respect to the metric (^). Thus, a;^ denotes 
the energy eigenvalue in units of Airf'T^/g. 



18 



We solved the system (P^, ( P^D for several values of the parameters rjQ, k, X and (?^$ and 
found only solutions with u"^ positive. For fixed values of the parameters we found several discrete 
normal modes which can be characterized by the number of nodes N of the fluctuation function 
Their eigenvalues a;^ increase with the number of nodes N. The lowest positive eigenvalue 
ujI corresponds to one node of the function see Fig. 8 (inlet). It seems to be likely that the 
system possesses an infinite number of discrete positive eigenvalues, forming a sequence 

with convergence to rj^. 





rio = 0.032 (A) 


r]o = 0.026 (i) 


N 


ujI (xlO-3) 


fit 




fit 


1 


0.9767 


0.9608 


0.66033 


0.65785 


2 


1.0098 


1.0080 


0.67165 


0.67144 


3 


1.0173 


1.0169 


0.67400 


0.67400 


4 


1.0201 


1.0200 


0.67486 


0.67486 


5 


1.0215 


1.0214 


0.67527 


0.67527 


6 


1.0222 


1.0222 


0.67549 


0.67549 


7 


1.0227 


1.0227 


0.67563 


0.67563 


8 


1.0230 


1.0230 


0.67571 


0.67571 



Table 1 

The eight lowest positive eigenvalues for rjQ = 0.032 and tjq = 0.026 on the sphaleron branch A 
for K = 0.4, A = and g,,^ = together with the fitted eigenvalues. 



For large N the eigenvalues uj^ can be well approximated by the formula 

^^=(^0-^)', (38) 

where C depends on the parameters rjo, k, A and (7^$. In Table 1 we give the first eight eigenvalues 
together with their fitted values for rjo = 0.032 and rjo = 0.026 on the branch A for fixed 
n = 0.4, A = and g^'p = 0. For the constants C we found C(?7o = 0.032) = 0.001 and 
C{rio = 0.026) = 0.00035. 

In Fig. 8 we show the lowest positive eigenvalue ujf as a function of rjo for k = 0.4, A = and 
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g^q> = 0. For all branches B, A and A the eigenvalue is a monotonically increasing function of 
?7o and vanishes in the limit rjQ ^ 0. 




5.2 The system {^^a^/^^/} 

For the system {'^a'^h'^f} the gauge zero modes are absent and we can calculate the differential 
equations and the static Hamiltonian directly. The system of differential equations is given by 



3a2 - 1 



+ /i2 + sinV + 4«'sin7(/'2 + 



3a2 sinV 





'2ah 




+ 






X 





S^^a/'sinV 



X 



+- 



X 



a cos/ sin/ + An'^a sin/ ( cos//'^ — 



sin//' 2a^ cos/ sin^/ ' 



+ 



x^ 



(39) 
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2a? 

— + gl^sinf + AK\'ih''-ril] 













X 





2gi^cosfsmfh f{AO) 



X 



—— (a cos/ sin/x*^ + 32/? '^a^ sin^/(a^ cos/ sin/ + a' f. 
xG ^ 

+4/t^ sin/x^ (2a^ cos/ sin^/ — a cos//'^x^ + 2a sinff'x — a' sin//' 
-5f^$acos/sinV/i^a;^))] 

cos/ sin fhx^ 



X 



+ 



x^G 



asm ff'x 
S/t^a^ sin^/ + 



K^a^ sin / + x^ 



(x^(2a2 + f^2^/iV)(l-2sin2/) 



+64fi:^a2 sin^ {a sin2/(a^(l - 2 sin^) - 2) + 2a cos/ sin//'x 
+a/'V + 2a'xsinV} 

+8K^ax^ |a sinV(a^(l - 4 sin^) - 2) + 6a cos/ sinff'x + 2a sin^ fx'^f 
—ax^f'"^ — 4 cos/ sin/a'/'x^ + 2 sin^/a'x — S'^.j.a sin^//i^x^|^ 



16^ 



, sin/ (^(x^ + 8fi;^a^ sin^/)(a(sin/ — cos//'x) — sin/a'a 



with 



G = x^ + IG/t^a^ sin/2(x2 + An'^a^ sin/^) . 
We now find for the dimensionless energy of the solutions of the differential equations 



^/ (41) 
(42) 



1 + 



8/t a sin / 



— - > dx . 
2 I 



^From this form we can define an appropriate metric for the fiuctuations \E' by 



1 + 



SK^a^ sin / 



x^ 



(ix 



(43) 



(44) 



Assuming the normalisiation of the fiuctuation functions according to (^Ij), we find the energy 
of the solutions to be Eahf = i^^, i- e. a;^ is again the energy eigenvalue. 

For the system {\l/a\l//j\E'/} we found normalizable solutions only on the saddlepoint branch 
A. For these normal modes the eigenvalue u"^ is negative and vanishes at the bifurcation points. 
Hence these normal modes represent an instability mode of the saddlepoint solutions. 

In Fig. 9 we show the negative eigenvalue uj"^ as a function of tjq for k, = 0.4, A = and 

S'tt* = 0. 



6 Summary and Discussion 

We have studied a model combining the Georgi-Glashow and the Skyrme systems interacting 
mainly through the so{3) gauge field, with an additional interaction term between the Higgs and 
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chiral fields. This interaction term, which breaks the chiral symmetry and fixes the asymptotic 
value of the chiral field at infinity, exploits the non vanishing VEV of the Higgs field in an 
essential way. 

The main emphasis of the study is the numerical analysis of the classical solutions of the 
model, which may be relevant to the semiclassical approach to monopole catalysis of Baryon 
decay . The most interesting feature of these solutions, both intrinsically and for the physical 
reason just mentioned, is the particular bifurcation patterns that they exhibit. These patterns 
are connected to similar bifurcations in the solutions of the gauged Skyrme model, not involving 
a Higgs field, which has led us to make a systematic study of the relation between the solutions 
of these two models. This has been carried out by considering particular limits, in terms of the 
independent parameters involved in the Monopole-Skyrme model, in which the solutions of this 
model merge with the solutions of the Higgs independent gauged Skyrme model. 

The above mentioned technical investigations form the centre of gravity of the present work. 
From the results obtained, some interesting observations of physical relevance can be made. 

One is the particular shape of the bifurcations the solutions exhibit, namely what we have 
referred to as 'butterfiy' in the text. These are reminiscent of the kind of bifurcations appearing 
in first order phase transitions. Unfortunately, in the present form of the model considered, we 
have not been able to make a concrete description for this phenomenon. This aspect of the work 
is under consideration. 

The other observation can be made more quantitatively. It was shown that in some limit 
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of the parameters, the Baryon resides in the core of the monopole while in the other hmit the 
converse, namely that the monopole resides inside the Baryon. In particular, for the values of the 
parameters fixed by the phenomenological values of the Higgs VEV and the pion decay constants, 
exactly one half of the Baryon charge resides inside the monopole core, and the rest outside. 
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